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ABSTRACT 

We report on observations of the Sunyaev-Zeldovich effect towards X-ray 
ROSAT clusters taken with a double channel (1.2 and 2 mm) photometer 
installed at the focus of the 15m SEST antenna in Chile. This paper describes 
the first results obtained for the high-z clusters S1077, A2744 and S295. 
Marginal detections were found for A2744 and at 1 mm for S1077. We discuss 
these data in terms of contamination of sources along the line of sight and give 
a constraint on the amplitude of the kinematic effect. 

Subject headings: cosmology: cosmic microwave background 
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1. Introduction 

The Sunyaev-Zeldovich (SZ) effect is a characteristic signature on the Cosmic Microwave 
Background (CMB) radiation spectrum originated from (inverse) Compton scattering of the 
photons by hot electrons in clusters of galaxies. There has been considerable interest in its 
detection because of its potential as a diagnostic tool for observational cosmology [physical 
process in the earlier universe; determination of Hq and qo and peculiar velocities of 
clusters) and for the study of the intracluster medium (see the original papers by ^eldovicE 



|fc Sunyaev 1969| , |Sunyaev fc Zeldovich 197!^ , [Sunyaev fc Zeldovich 19811 ). Its detection has 



been the goal of many observational searches (e.g. Birkinshaw 1990|) , mostly in the Rayleigh 



- Jeans (R-J) part of the spectrum, where the scattering leads to an intensity decrement. 
However, interpretation of measurements in the R-J side is plagued both by the possible 
radio emission from sources within the clusters and by the smallness of the effect. After 
many attempts to detect this effect (see e.g., the review paper by Birkinshaw 1990) radio 
observations show the expected decrement at centimeter wavelengths towards A2218, A665 
0016+16, A773 and Coma ([Birkinshaw 1991| , [Klein et al. TMl] , |J ones et al. 1993| , [Grainge| 



et al. 1993[ , [Herbig et al. 1995[) and at 2.2mm towards A2163 ( [Wilbanks et al. 1994| ). 



Measurements near the planckian peak and on the Wien side are in principle more attractive 
since: (a) the intensity enhancement relative to the planckian value is larger than the 
magnitude of the R-J decrement; (b) the simultaneous detection of the enhancement 
(positive) and decrement (negative) on the same cluster provides an unambiguous signature 
of their presence; (c) sources in the cluster are expected to give a negligible contribution at 
high frequency; (d) because of the large bandwidth the sensitivity of bolometer systems is 
excellent. 

We have built a photometer with two bolometers centered at 1.2 and 2 mm to feed the 
O.A.S.I. (Osservatorio Antartico Submillimetrico Infrarosso) telescope installed at the 
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Italian base in Antarctica (|Dall'Oglio et al. 1992 ). This photometer was adapted to the 
focus of the S.E.S.T. and its performances were tested during an observing run in September 
1994. Most of the allocated time was spent in installing the photometer and testing the 
equipment (pizzo et al. 1995a|) but some useful observations were gathered towards S1077, 
A2744 and S295. 



The Instrument 



The instrument was described in detail elsewhere ( |Pizzo et al. 1995a|) , here only the 



main characteristics are summarized. The photometer uses two Si-bolometers cooled to 0.3 
K by means of a single stage ^He refrigerator (|Pizzo et al. 1995bD . Radiation coming from 



the telescope is focussed by a PTFE lens and enters the cryostat through a polyethylene 
vacuum window and two cooled filters. A dichroic mirror (beam-splitter) at 4.2 K divides 
the incoming radiation between two f/4.3 Winston Cones cooled at 0.3 K located orthogonal 
each other. The wavelength ranges are defined by two interference filters centered at 1.2 
and 2 mm cooled at 4.2K with bandwidth 350 and 560 /im respectively. The 2mm band 
includes the peak brightness of the decrement in the S-Z thermal effect, while the 1.2mm 
bandwidth is a compromise between the maximum value of the enhancement in the S-Z 
and the atmospheric transmission. 

Since our goal is the simultaneous observation of the same source region, comparable beam 
sizes at the two different frequencies were achieved using Winston cones with the same 
entrance sizes for both channels, selected to have roughly a diffraction limited configuration 
for the 2mm channel. Beam shape and dimension were measured using raster scans on 
Saturn, Jupiter and Uranus (details on the procedure used and on the map reconstruction 
are in Haslam (1974), Sievers et al. (1991), Lemke (1992)). The planet maps show a good 
alignment (within 2") of the two beams and their symmetry around the optical axis. Their 
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Half Power Beam Widths turn out to be 44" (1.2 mm channel) and 46" (2 mm channel). 

Since clusters of galaxies are extended sources, in order to have the reference beam as far 
as possible from the peak region of the X-ray emission, the chop throw was set to the 
maximum possible amplitude, i.e. that corresponding to a beam separation on the sky of 
135". The coupling between the antenna and the photometer was checked by comparing 
the output signals from known sources at the 1.2mm channel of our photometer with those 
at the 1.2mm MPI bolometer available at SEST ( [Kreysa 1990|) , both having a filter very 



similar in transmission and shape. Taking into account the correction factors due to the 
different fields of view, equal signals within the errorbars have been found. 

The detector noise of both channels measured at the focus was about h^nV / \/Hz. 
Responsivities were inferred with Venus, Saturn, Uranus, Mars and Jupiter taking the 
planet brightness temperatures from Ulich (1981 and 1984): the effective responsivities at 
the focus turned out to be 1.8/iV/K and 0.7/iV/K at 1.2 and 2 mm respectively. This means 
that the measured values for the N.E.T.s were 28 and 70 mK/V Hz, i.e. 14 and 35 mK 
in Is integration time. 10 % is the average uncertainty on these figures. Converting these 
values in thermodynamic temperatures gives for the measured sensitivities ^ = 0.014 and 
= 0.02 at 1.2 and 2 mm respectively (Is integration time). 

Atmospheric transmission was monitored by frequent skydips, i.e. telescope scans at 
constant azimuth and different elevation are carried out to compare the atmospheric 
emission with that of an absorber moved into the path of the reference beam. 

The position of the sources was found by first pointing a nearby radio-loud quasar with 
strong millimetric fluxes and then the X-ray coordinates. With this procedure we reduce 
the pointing errors to less than 2". 



- 7- 



3. The observing strategy 

Targets were selected for their high X-ray luminosities in the ROSAT band (O.l-i-2.4 
keV). They are extended X-ray sources identified as optical high-redshift {z ~ 0.3) clusters. 
Names, redshifts and ROSAT luminosities are listed in table 1. If we assume a core radius 
of 250-i-400 kpc in standard cosmologies their apparent size is 45-i-70", a chop throw of 
135" means that at the reference beam position the ratio between the electron density, 
ne{e = 135"), to its central value, 71^(6 = 0), is ~ 0.28 ^ 0.42, i.e. we lose 30^40 % of the 
signal because of the limited chop throw. The effective values of y, taking account of the 
SZ effect in the off-center positions, can be roughly estimated only for A2744 for which a 
deep PSPC exposure is already available[^ The gas density was derived from the best-fit 
value of a /5-model to the surface density profile and the temperature was assumed to be 7 
keV. y turns out to be 1.3 10""^. For S1077 and S295 the X-ray properties listed in table 1 
were taken from the ROSAT All Sky Survey and the luminosities were estimated assuming 
Te = 7 keV. 

We have checked that the selected candidates do not contain IR sources of the 
IRAS Faint Source Catalogue ( pVIoshir et al. 1989|) . But sources fainter than the 90 % 
completeness limit of this survey (~ 230 mJy at 60 fim) may fall in our beam. On the basis 
of our observations on galaxies performed at SEST and IRAM (Andreani fc I'ranceschinl 



1992| and [I'ranceschini fc Andreani 1995| ), the expected millimetric emission, scaled at the 



cluster redshift, is less than 0.8 mJy at 1.2 mm. This fiux would produce an output signal 
in antenna temperature of 0.03 mK, which is lower than the expected signals from the S-Z 
effect (see below). 

The magnitude of the expected signals from the effect was estimated by convolving 



^"ROSAT HRI observations are already scheduled for S1077 and S295 
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— y {x — 4) with the spectral filter response, the atmospheric transmission and the 
beam-shape, where T is the CMB temperature, x = hu/kT and y = j'{kTe/mc^) Ueardi 
is the comptonization parameter, ng, Tg being the electron density and temperature. We 
find: ^ ~ 0.69y, ^ ~ — 1.77y at 1.2 and 2 mm respectively for the thermodynamic 
temperatures. 

In order to reduce the atmospheric noise and systematics from the antenna two 
combined observing strategies were exploited. 

(1) To reduce temporal and spatial drifts of the atmosphere a three beams technique (beam 
switching + nodding) was applied. Let us call position A that with the source in the right 
beam and the reference beam on the left and position B that with the source in the left 
beam with the reference beam on the right, the antenna tracks each position to integrate 10s 
and moves according to the following temporal sequence ABBAABBAA....A. 20 positions 
are tracked so that the entire cycle lasts 200s, which is the shortest integration time and 
which in the following will be called one scan (see §4 for further details). The duty cycle 
was 80 %, i.e. additional 4s are spent for each movement of the antenna between position 
A and B. 

(2) In order to properly quantify the contribution to the source signals from unknown 
systematics we measured them using the following procedure. Each source was integrated 
over time chunks of 600s (3 of the above scans), this time interval plus the needed overheads 
gives a total tracking time on the source of 15 minutes. The same time was spent on a 
blank sky located in equatorial coordinates 15 minutes ahead in right ascension. This 
means that the antenna tracks twice the same sky position in horizontal coordinates: once 
ON the source and the other OFF the source. This enables us to compare the two different 
measurements and quantify the spurious signals (atmosphere + antenna systematics). 
The choice of 600s of integration ON and OFF the source is a compromise between the 
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minimization of the time wasted on overheads and the need of minimizing the atmospheric 
variations between one observation and the other. 



4. Data Analysis 

As described in §3 the data consist of chunks of 3 scans, each containing 20 sub-scans 
of 10s corresponding to either position A or B. The reduction procedure is done in several 
steps. (1) One scan data (200s: 440 data points) are roughly Gaussian distributed with 
a width larger than the detector noise because of the atmospheric noise. In order to get 
rid of individual data inconsistent with the underlying Gaussian distribution, each datum 
with a value exceeding 3 times the r.m.s. of the entire scan was rejected. This criterion 
eliminates spikes due to equipment malfunction or atmospheric sudden variation. The 
fraction of rejected data is 1%. (2) Within each sub-scan high frequency atmospheric noise 
was filtered with the Savitzky-Golay algorithm (Press and Teukolsky, 1990), which smooths 
the time series within a moving window by approximating the underlying function by a 
polynomial of 2"'^ 3^*^ order. At each point the least-squares fit to the n^ + ur + 1 
points in the moving window is assigned, where n^ and n^^ are the parameters chosen to 
smooth the data over a bin of Stg {to being the integration time of the lock-in amplifiers, 
i.e. twice the inverse of the chopping frequency, 4.4Hz). A mean value of the differential 
sky temperature (T^ or Tg) is then obtained for each sub-scan, whose variance is found 
with a procedure of bootstrap resampling (e.g.. Barrow et al. 1984). This method is widely 
used to take into account correlation among data (induced in this case by the filtering). (3) 
The signal is obtained by subtracting each couple of subscans AT^ — ^ with variance 
2af = a\ + a%. (4) For each 200s scan the average value is AT^ = * , with a 

scatter of - E,:(ATf-AT;^)^a-^ 
scatter or - (^.a-^)(M-i) ■ 
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The weighted mean over the 3 scans is (at the sky position P) ATp = — "'.g'" , with a 
scatter about the mean of Up = -^"=^'3^ — zf^ — ^. 



For each AT^ a weight is assigned: = (a^ + dp) ^, where represents the variance 
due to short-term fluctuations, while Up that due to medium-term variations. 

The final weighted mean, over the sky position P, is given by: 



y-3 ^rjpa 



y 



with estimated variance: 



,_i E^=i(Ar°-(Ar%f -^^ 

Values obtained when the antenna tracked the source for 600s (3 scans), hereafter called 
signal ON, are subtracted from those taken when the antenna was 15 minutes ahead in 
R.A. (hereafter called signal OFF): 



l^T^sz = (AT;)oAr - (at;)opp 

and the quadratic sum of the two standard deviation has been used as errorbar: 

^"sz — {(^'f)oN + (c/)off- {^T^)on and {ATf)oFF with their errorbars are plotted in figure 

1 for the three clusters and for each channel. 

The final values, the weighted means over all the subtracted scans in thermodynamic values 
and the corresponding y parameters are listed in table 1. 
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5. 



Discussion 



The size of the y value of A2744 derived in table 1 disagrees with the estimated one 
(see §3) and cannot be reliably ascribed to the sole S-Z thermal effect unless some relevant 
information is missed in the X-ray data. The resulting data towards A2744 at 1.2 and 2 
mm and S1077 at 1.2 mm could be therefore originated from other sources. Here possible 
contaminations are briefly discussed. 



Peculiar velocities of clusters could contribute significantly to the signals for values 
larger than Vr > 1000 km/s. In fact the ratio between the S-Z kinematic effect, ^ = ' t 
(where t — J arndl is the optical depth for Thomson scattering along the hne of sight and 
Vr is the peculiar velocity), and the thermal effect depends on the cluster peculiar velocity 
and the gas temperature. The lack of any other observations of peculiar velocities on these 
clusters makes the quantification of this contribution impossible and on the basis of the 
present knowledge we cannot exclude that a fraction of the signals is due to it, even if large 
peculiar velocities seem quite unlikely in most cosmological scenarios. The present data 



can be used to constrain the combination r • — which turns out to be < 3 10 ^ for A2744, 

c ' 

S1077 and S295 (3(7 values). 



As mentioned in §3 a strong contamination from sources in the cluster seems to be 
unlikely: (a) clusters do not contain many spirals, (b) late-type galaxies at this redshift 
would have a millimetric emission below than the expected signals from the S-Z effect and 
(c) we are looking at the very center of the cluster which, as shown from optical images. 



5.1. 



Contamination of kinematic effect 



5.2. 



Contamination of cluster sources 
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is devoid of sources. Only the cluster S295 has three ellipticals very close to the center of 
the X-ray image and part of their emission can therefore falls in our beam. If there was a 
significant non-thermal emission we would expect to find positive signals in both channels, 
contrary to our results (see table 1). 

Let us suppose that one or more sources fall in one of the reference beams. In this case if 
a mm-emitting source was located in the right beam, we would find a positive signals in 
both channels. If it was in the left reference beam the signals would be negative in both 
channels. This seems to be unlikely since we do not find these systematics (see table 1). 
However contribution from cluster sources deserves further investigations and will be the 
goal of future projects. 

It is also unlikely that both signals are contaminated from the emission of unresolved 
sources: the estimate of their confusion limit in the beam made by Franceschini et al. 
(1991) gives at 1.2 mm AI ~ 0.4 mJy, i.e. AT ~ 0.06 mK. 



5.3. Contamination by other spurious signals 

We have considered other effects affecting the signals: (a) the selected clusters have 
sizes comparable to our beams and beam dilution cannot significantly lower the detection 
probability ( [Rephaeli 19871 ), but fiuxes are reduced up to 30-i-40 % because of the limited 



chop throw, (b) Diffraction patterns could enter the beams and alter the signals but 
constant systematics of this type should be removed by our observing strategy (see §3). 
However, if the antenna did not track precisely the same paths relative to the local 
environment, because of a loss of synchronization, some spurious signals survive and we 
cannot exclude at present that this never happened, (c) We checked if sources are located 
in the OFF position (15min ahead in R.A.): from the IRAS Faint Source Catalogue and 
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5GHz NRAO survey it turns out that none fall in these areas. The IRAS 100 jim sky maps 
do not show any emission from galaxy cirrus at a level greater than their zero level. 

The authors are indebted to Glenn Pcrsson and to the ESO Workshop team at La Silla. 
This work has been partially supported by the P.N.R.A. (Programma Nazionale di Ricerche 
in Antartide) and made use of the Skyview Database, developed under NASA ADP Grant 
NAS5-32068. Very useful comments by an unknown referee helped to improve this paper. 
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ROSAT Cluster 

name redshift Ljs:(0.1^2.4kcV) 
(1045 erg/s) 


i^T) therm 

1 mm 2 mm 
(mK) (mK) 


Compt. par. 

(10-^) 


S1077 0.312 (0.66d=0.16) 
A2744 0.308 (2.00±0.05) 
S295 0.299 (1.30±0.15) 


(0.9 ± 0.3) (-0.8 ± 0.5) 
(0.8 ± 0.3) (-1.4 ± 0.4) 
(0.1 ± 0.3) (-2.5 ± 1.0) 


(5.3±1.8) 
(3.8±1.2) 
< 5.3 (3(7) 
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Fig. 1. — Antenna temperature differences for positions OFF (empty squares) and ON (filled 
squares) the source. Right panels refer to 2 mm channel, while the left ones to the 1 mm 
channel. Note that in most cases the 2 mm signals ON the sources lie below the corresponding 
one OFF the sources. The statistical significance is however very low. The corresponding 
panels for the 1 mm channel show the opposite trend meaning that the signals from the 
sources are higher than those OFF the sources. 
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